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Abstract 

We calculate the solar-modulated energy spectra of cosmic-ray antiprotons (p's) 
from two candidate primary sources, i.e., evaporating primordial black holes and 
the annihilation of neutralino dark matter, as well as for the secondary p's pro- 
duced by cosmic-ray interactions with interstellar gas. A large enhancement toward 
the solar minimum phase emerges in the low-energy flux of p's from the primary 
sources, whereas the flux of the secondary p's, falling steeply below 2 GeV, does not 
significantly vary. This enables us to conduct a very sensitive search for primary 
p components by precisely measuring the p spectrum, especially at low energies, 
throughout the forthcoming solar minimum phase. 



1 E-mail address 

2 E-mail address 
3 E-mail address 



mitsui@icepp.s.u-tokyo. ac.jp 
maki@icepp . s . u-tokyo .ac.jp 
oriton@icepp.s.u-tokyo. ac.jp 



1 Introduction 

Since the discovery of antiprotons (p's) in cosmic rays fll], ^ [| , their origin has attracted 
much interest (for a review, see Ref. 0). From a standpoint of cosmology and elementary 
particle physics, it would be of great importance to detect or to set a stringent limit on 
the primary p components from such sources as the annihilation of neutralino dark matter 
P, [], |7], 0], evaporating primordial black holes (PBHs) || [TD| and superconducting cosmic 



strings ||11|| . No evidence for such primary p components has been obtained so far. 



For the very sensitive search for primary p's anticipated in the near future experiments 



I4 , a potential source of background could be the secondary p's produced by cosmic-ray 



interactions with interstellar gas, since the models of the cosmic-ray propagation predict 



the secondary p flux at a level close to the present observations j|, 13, 14, 15, 16|. However, 
the energy spectra of the primary p's are generally expected to be rather soft (i.e., flux 
increases with decreasing energy), whereas the spectra of the secondary p's should steeply 
fall below 2 GeV due to the kinematics of the p production. The signals of the primary 
p's can, therefore, be searched for in the low-energy regions. 

In this Letter, we point out that such low-energy portion of the primary p's should show 
a large enhancement toward the solar minimum phase and a rapid decrease afterward, 
whereas the secondary component, falling steeply below 2 GeV, does not significantly 
vary. The next solar minimum is expected to arrive in a few years: From the observations 



of 22-yr solar cycle over the past 45 years by, e.g., the CLIMAX neutron monitor fl7|| , 
one can expect that the next solar minimum should be arriving in the years 1997 to 
1998, which should be followed by a rapid increase of solar activity. The low-energy flux 
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of the primary p's, e.g., those from evaporating PBHs, would increase by a factor of 5 
in 1997-1998 as compared to the level in 1993, and would then rapidly decrease by a 
factor of 15 in 2000. On the other hand, the secondary p component shall remain stable 
within a factor of 1.5 throughout this period. Therefore the extraction of the primary 
p components would become possible with precision measurements of the p spectra over 
the next several years, providing a unique opportunity either to detect or put a stringent 
limit on the primary p components. 

2 Evaporating primordial black holes 

Primordial black holes (PBHs) may exist as a result of initial density fluctuations, phase 
transitions, or the collapse of cosmic strings in the early Universe (for a review, see 
Ref. fL8|). Black holes emit particles and evaporate by quantum effects (Hawking radia- 
tion) jn|; in particular, PBHs are the only ones with a mass small enough for the quantum 
emission rate to be significant, possibly producing an observable flux of p's [[| |It]]. 

In Ref. [[TIJ, we investigated low-energy cosmic-ray p's from evaporating PBHs using 
the Monte Carlo simulation code JETSET 7.4 to obtain the source spectrum of p's 
from PBHs. The local interstellar flux of these p's was then calculated using a 3-D Monte 
Carlo simulation based on the diffusion model of the cosmic-ray propagation. The spatial 
distribution of PBHs was assumed to be proportional to the isothermal distribution of 
halo dark mass, i.e., ph{r) = Ph©(^0 + r c)/( r2 + r c) HHL where r is the Galactocentric 
distance with r = 8.5 kpc being the position of the Solar system, r c = 7.8 kpc is the core 
radius, and ph© = 0.3 GeV cm -3 is the local halo density. Our 3-D simulation showed 
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that only PBHs within a few kpc of the Solar system contribute substantially to the local 
interstellar p flux (see also Fig. [I]). We then found that the local PBH explosion rate, 1Z, 
must be less than 1.7 x 1CU 2 pc~ 3 yr _1 in order not to conflict with the BESS '93 data ||. 
Nevertheless, we will show later that future experiments can have sensitivity to TZ value 
down to ~ 1/30 of the above upper limit. In this Letter, we use the local interstellar p 
flux obtained in Ref. [|nj after scaling it for an appropriate value of TZ. 



3 Annihilation of neutralino dark matter 

Supersymmetric (SUSY) models for particle physics have been extensively studied in 



recent years |22fl . In most models, the lightest neutralino x 1S taken to be stable and 
thus serve as cold dark matter (CDM) in the Universe. If the neutralinos comprise the 
Galactic halo, their annihilation may produce a detectable flux of p's || ||, [7], ||, since a 
pair of neutralinos annihilates into quarks, leptons, gauge bosons or Higgs bosons, which 
then decay or fragment into various particles including p's. 

Here we study this process in the framework of minimal N — 1 supergravity models 
with radiative breaking of the electroweak gauge symmetry In this framework, one 



can calculate all SUSY particle masses and couplings by solving the renormalization group 
equations with only four SUSY soft breaking parameters: m , m 1 / 2 , tan/3 and A . The 
higgsino mass parameter p is also determined, except for its sign, from the requirement 
of the electroweak gauge symmetry breaking. We assume p to be positive, and the top 
quark mass to be m t = 175 GeV. Part of the parameter space has already been excluded 
by theoretical and experimental bounds, including recent results of searches for SUSY 



particles at Tevatron [ 24L EH] and LEP 1.5 [26 



We then compute the value of Q x h 2 , where Q x is the neutralino relic density in 
units of the critical density of the Universe and h is the Hubble constant in units of 
100 km s _1 Mpc _1 . Recently, Berezinsky et al. [^71] have pointed out that the expected 
value of Q x h 2 falls within the range of 0.2 ±0.1 in most cosmological models from the 
viewpoint of the age of the Universe. Thus, we choose here four representative parameter 
sets which give Q x h 2 = 0.18 (see Table [I]). Neutralino dark matter with this value of 
Q x h 2 is most likely pure bino, and annihilates predominantly into bottom quarks (bb) or 
tau leptons (r + r~). 

To obtain the source spectrum of p's from the neutralino annihilation, we use the frag- 
mentation functions extracted from JETSET 7.4 p0 |. Since most of the p's are produced 



as the fragments of bb (no p from t + t~), the shape of their source spectrum does not 
depend much on the parameter set, although its absolute value scales with SB b i, where 
S is the neutralino annihilation rate per unit volume and B b i is the branching ratio into 
bb. The local interstellar flux of these p's is then calculated using a 3-D Monte Carlo 
simulation based on the diffusion model, as in the case of evaporating PBHs. Note that, 
whereas the emission rate from evaporating PBHs per unit volume is proportional to their 
density, S is proportional to the squared neutralino density; thus the local flux of p's from 
the neutralino annihilation is rather sensitive to the halo density distribution. 

We will show later that the flux of p's from the neutralino annihilation is usually too 
small to be observed under conventional astrophysical assumptions, i.e., if the neutralino 
has the homogeneous, spherical and isothermal distribution with p^o = 0.3 GeV cm -3 . 
There are, however, at least three astrophysical uncertainties which possibly increase the 
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p flux by a factor of £ ^> 1: First, the Galactic halo may be flattened toward the Galactic 



disk | 28fl . This shall increase the local halo density /?h© by a factor of ~ 2 over that of 



a spherical halo ||29|| , resulting in the enhanced annihilation rate by a factor of £ ~ 4. 
Secondly, a non-dissipative gravitational singularity (NGS) may reside at the Galactic 
center (GC) P0| , resulting in the halo density distribution being ph(r) = ph©(?"/ ? ''©)~ 1 ' 8 
for r > 0.1 pc. If this is the case, the neutralino annihilation rate would be very high at the 
GC. Figure [I] shows the r-distribution of the primary sources (the neutralino annihilation 
and evaporating PBHs) which contribute to the integrated p flux at the Solar system. It 
can be seen that, whereas only the annihilation occurring within a few kpc of the Solar 
system would contribute in the case of the isothermal density distribution, in the NGS 
model, p's from the annihilation at the GC would dominate the local interstellar p flux, 
leading to £ ~ 200. Thirdly, high-density CDM clumps with a density enhancement factor 
of 10 2 -10 9 might be generated in the early Universe [31j]. If a small portion (a few %) 



of neutralino dark matter is in the form of such clumps, the expected p flux could be 
enhanced by a factor of £ > 20 P2"| . Furthermore, if there exists such a clump in the 
vicinity of the Solar system, the resultant p flux at the Earth could be further enhanced 
by a large factor. 

4 Secondary p flux 

When high-energy cosmic rays interact with the interstellar gas, secondary particles are 
produced including p's. The expected flux of such secondary p's was calculated by a 



number of authors 0, 13, 14, 15] who showed that the expected energy spectrum of the 



secondary p's drops off steeply below 2 GeV. This feature originates from the kinematics 
of the p production, and is independent of the details of the cosmic-ray propagation. 



Recently, one of us has calculated the secondary p fluxes |16| using the two models 
of the cosmic-ray propagation, i.e., the Standard Leaky Box (SLB) and Diffusive 
Reacceleration (DR) models, by taking the most plausible input data such as the 
parent proton flux Q The resultant spectra shown in Fig. |2| have very similar shapes, 
although the absolute flux in the SLB model is 2.5 times larger than that in the DR 
model. 



5 Solar modulation 

When antiprotons go into the Solar system, their energy spectra are changed due to the 
diffusion, convection, and deceleration processes by the solar wind and the interplanetary 
magnetic field (solar modulation). The minimal model to describe the processes is the 
spherically symmetric model P5L [57|, in which all quantities are dependent only on the 



Heliocentric distance g. At the boundary of the Heliosphere (g = gb) the cosmic-ray 
energy spectrum is assumed to match the interstellar spectrum. The energy spectrum at 
the Earth (g = gE = 1 AU) can then be calculated by solving the diffusion-convection 
equation |37| using the numerical technique developed by Fisk 



As usual, we take the solar wind speed V = 400 km/s and the radius of the Heliosphere 



gb = 60 AU. Since it was shown |T6| that the spectrum at the Earth is not significantly 



affected by the position dependence of the diffusion coefficient k, we take k to be position 



1 The interstellar proton flux is taken to be 1.5 x 10 4 /3 1 P 2 74 m 2 sr 1 s 1 GeV 1 , where P is the 
momentum in units of GeV, based on the recent observation of the LEAP experiment [p5[ . 



independent and, as normally done, to be proportional to the particle rigidity R as well 
as to the speed (3, i.e., k = (3k\R, where the coefficient k% varies with solar activity. The 
solar activity level is then represented by a single parameter: (ftp = (g b — Qe)V/3ki, which 
roughly corresponds to the "0 parameter" in the force field approximation j37| Q, i.e., 
the energy loss in the Heliosphere per particle charge. With <pp varying between ~ 350 
and ~ 1500 MV from the solar minimum to the maximum, this model reproduces well 
35, [IB] the large variation (a factor of 20) of low-energy proton and helium fluxes 



as well as each energy spectrum at various stages of solar activity. Since the soft energy 
spectra of the primary p's are similar to those of proton and helium, we use this model 
of solar modulation to calculate the solar-modulated spectra of the primary p's. 

To forecast the time variation of (j)p, we can utilize the data of the CLIMAX neutron 



monitor |T7|], because the <pp parameter shows the excellent correlation with the count 



rate of the neutron monitor [p9(1 . The CLIMAX data over the past 45 years clearly show 
22-yr solar cycle. Especially the recent trend of the data shows strong similarity to those 
in ~ 1950. From this, one would expect the next solar minimum in the years 1997 to 1998 
with the 4>f parameter reaching ~ 350 MV, which should be followed by a rapid increase 
of solar activity (<j> F ~ 1000 MV) in 2000. 



6 Results and discussion 

Figure § (a) shows examples of the primary p fluxes as well as the two secondary p fluxes 
(SLB and DR) solar-modulated with <fip — 350, 550, and 1000 MV, which would corre- 
2 Analyses of proton data show that (ftp ~ <\> + 50 MV. 
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spond to 1997, 1995, and 2000 respectively. The uppermost curves (<pp = 0) represent the 
interstellar spectra. We have taken TZ = 5 x 10 -3 pc~ 3 yr _1 for the p flux from evaporating 
PBHs, and the case #1 with £ = 25 for the p flux from the neutralino annihilation. It 
should be noted that the secondary p fluxes are rather insensitive to the solar mod- 
ulation. This is because the interstellar fluxes decrease with decreasing energy below 2 
GeV so that the two effects of the solar modulation, i.e., deceleration and flux suppression 
cancel each other. Contrastingly, the primary p fluxes would show a large enhancement 
toward the solar minimum in a very similar way to the proton flux. Note that even 
the "enhanced" solar minimum fluxes of the primary p's are very much suppressed as 
compared to the interstellar fluxes. 

Based on the present understanding of the cosmic-ray propagation, we expect that the 
secondary antiprotons will most likely exist. Therefore, we show in Fig. |2| (b) the solar- 
modulated (4>f = 350 and 1000 MV) p fluxes in three cases, i.e., (i) secondary (in the SLB 
model) only, (ii) secondary (SLB) plus p's from evaporating PBHs, and (iii) secondary 
(SLB) plus p's from the neutralino annihilation. Figure ^| (c) shows the corresponding 
spectra with the DR model of the secondary p's. 

The unknown absolute intensities of the primary p fluxes (TZ and £ for PBHs and 
neutralino dark matter respectively) were chosen here not to conflict with the existing 
observations. Note that the value of TZ taken here (see above) is 1/3 of the upper limit 



obtained in our previous paper |T0| using the BESS '93 data ||. Such a weak signal, how- 
ever, will appear as a striking enhancement of the low-energy p flux at the solar minimum. 
Also shown is the expected statistical accuracy estimated for a future measurement by the 



BESS experiment with a total 48-hour flight [12|]. As seen in these figures, the existence of 
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novel primary p sources at this level could be detected without any ambiguity It can also 
be seen that the signals of the primary p's will rapidly diminish after the solar minimum 
with the increase of solar activity. This strong variation with solar activity would provide 
further evidence for the primary p's, because the spectra of secondary p's are not largely 
affected by the solar modulation. 

In the past, the observational data as well as theoretical calculations were usually 
presented in terms of the p/p ratio. However, as shown in Fig. [| the p/p ratio is rather 
insensitive to the solar modulation of the primary p's, since the denominator of proton flux 
varies in a similar way to the primary p components. Therefore, one has to measure the 
absolute p flux and its spectrum in order to fully utilize the characteristic solar modulation 
of the primary p's. 

Finally we discuss the attainable sensitivity to the primary p components expected in 
the future experiments. Figure [| (a) and (b) show the p spectra corresponding to Fig. 
(b) and (c), respectively, for much smaller contributions of the primary p's. As seen in 
the figures, these small signals can also be detected with the expected accuracy. The 
corresponding sensitivity to the PBH explosion rate 72. is 1 x 10~ 3 pc~ 3 yr _1 in case of the 
SLB model, and 5 x pc _3 yr _1 in case of the DR model. These values are respectively 
1/17 and 1/34 of the best existing upper limit derived from the BESS '93 data. Note that 
the explosion rate 1Z = 1 x 10~ 3 pc~ 3 yr _1 is 9 orders of magnitude lower than the upper 



limit obtained from searches for TeV 7-ray bursts [41] , and is a few hundreds times lower 



than the upper limit derived from the anisotropy of the diffuse 7-ray flux [42]. On the 



standard assumption of the PBH initial mass spectrum [1C], this value of 1Z corresponds 
to the average PBH density in the Universe of flpBn = 3 x 1CT 10 , which is a factor of ~ 30 
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below the upper limit deduced from the cosmological diffuse 7-ray background [f43| . 

Figure ^ also shows that, in order for the p flux from the neutralino annihilation to 
be observable for the case #1, the enhancement factor £ must be greater than 8 and 4 
when the SLB and DR models, respectively, are used for the secondary p's. Since such 
enhancement is possible within astrophysical uncertainties as discussed in Section 3, p's 
from the neutralino annihilation could be observed in the near future. 

In conclusion, if future experiments allow us to precisely measure the low-energy 
cosmic-ray p flux at the solar minimum, we should be able to conduct a very sensitive 
search for primary p components by utilizing the characteristic solar modulation of the 
primary p's. 
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Figure captions 



Figure 1: Distribution (along the Galactocentric distance r) of the novel primary sources 
(the neutralino annihilation and evaporating PBHs) contributing to the integrated p flux 
F at the Solar system, which is located at r = 8.5 kpc. 

Figure 2: (a) The expected energy spectra of p's from evaporating PBHs with 1Z = 
5 x 10~ 3 pc~ 3 yr -1 (dashed lines) and from the neutralino annihilation for the case #1 
with £ = 25 (dotted lines), as well as the secondary p's in the SLB (thick solid lines) and 
DR (thin solid lines) models. The curves correspond, from top to bottom, to 4>f values of 
(interstellar), 350, 550 and 1000 MV. (b) The expected spectra for the secondary p (SLB) 
only (solid lines), secondary p (SLB) plus p's from evaporating PBHs (dashed lines), 
and secondary p (SLB) plus p's from neutralino annihilation (dotted lines). The upper 
and lower curves correspond to <f)p = 350 and 1000 MV respectively. The normalization 
parameters for the primary sources are the same as in (a). Also shown is the expected 



statistical accuracy of a future observation [12]. (c) Same as (b) with the DR model of 
the secondary p's. 
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Figure 3: The p/p ratio corresponding to Fig. § (b) and (c). The lower and upper curves 
correspond to <pF = 350 and 1000 MV respectively. The proton spectra are obtained by 
solar-modulating the interstellar spectrum: 1.5 x 10 A /3~ 1 P~ 2 ' 74 m~ 2 sr~ 1 s~ 1 GeV~ 1 f3~5fl , 
where P is the momentum in units of GeV. 

Figure 4: (a) Same as Fig. |2] (b) but for TZ = 1 x 10~ 3 pc _3 yr _1 and £ = 8. (b) Same as 
Fig. @ (c) but for TZ = 5 x 10" 4 p^yr -1 and £ = 4. 
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Tables 



Table 1: Representative sets of SUSY soft breaking parameters (assumed A = 0) and 
relevant quantities derived from them, i.e., (i) the neutralino mass, m x ; (ii) the mass of 
the lightest charged SUSY particle (LCSP), m LCS p, together with its identity (chargino 
xf or stau fi); (iii) the neutralino annihilation rate per unit volume in the local region 
of the Galactic halo, S = (pho/' m x) 2 ( cn; )h? where (<Jv)h is the thermally-averaged cross 
section for the neutralino annihilation in the Galactic halo; and (iv) the branching ratio of 
the neutralino annihilation into bottom quarks, B b i. Each set gives Q x h 2 = 0.18. Masses 
m o, m i/2, m x an d ^lcsp are shown in units of GeV, S in units of 10~ 32 cm^V 1 , and B b i 
in units of %. 
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